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Abstract. Operators in quantum mechanics - either observables, density or evolution 
operators, unitary or not - can be represented by c-numbers in operator bases. The position 
and momentum bases are in one to one correspondence with lagrangian planes in double 
phase space, but this is also true for the well known Wigner-Weyl correspondence based 
on translation and reflection operators. These phase space methods are here extended to the 
representation of superoperators. We show that the Choi-Jamiolkowsky isomorphism between 
the dynamical matrix and the linear action of the superoperator constitutes a ’’double” Wigner 
or chord transform when represented in double phase space. As a byproduct several previously 
unknown integral relationships between products of Wigner and chord distributions for pure 
states are derived. 


PACS numbers; 03.65.-w, 03.65.Sq, 03.65Yz 



Representation of superoperators in double phase space 

1. Introduction 


2 


The action of a unitary operator on states in quantum Hilbert space corresponds 
semiclassically to a classical transformation, x_ x+, in the corresponding classical phase 
space with points x = {q\, ...,qi^,pi, ...,Pn)- Alternatively, one may place the points x+ 
and x_ each in its own phase space so that the graph of a canonical transformation becomes 
a 2N- dimensional (2A-D) surface within the double phase space, built from the product 
space R^^ x R^^, with points X = (x+, x_) [[T]|. 

In strict analogy to the semiclassical correspondence of an (integrable) quantum state 
with a classical N-D surface in R^^ [|2l|3l, evolution of the operator corresponds classically to 
movement of the surface in its phase space. Thus, a coordinate transformation in phase space, 
x_ x+, implies a mapping X X' , in double phase space (for instance, a normal form 
transformation [|4l|3||) that may be seen to propagate the classical surface for a given canonical 
transformation in double phase space, corresponding to the action of a superoperator as it 
evolves ordinary operators in quantum mechanics. 

One should keep in mind some essential differences between the classical and the 
quantum scenarios. Above all there is the uncertainty principle: while vectors in Hilbert 
space correspond semiclassically to (complex) functions on the lagrangian manifolds on 
single phase space - say the position or the momentum basis wave functions- operators are 
represented by their matrix elements from a pair of such manifolds. When viewed in double 
phase space these pairs are themselves 2A-D lagrangian surfaces in 4A-D double phase space. 

However, some special coordinate transformations are allowed that take these obvious 
lagrangian coordinate planes of double phase space, that is, initial and final positions or 
momenta, to new coordinate planes that cannot be so decomposed. This is just the case 
of phase space labels underlying the Weyl representation (i.e. the Wigner function in the 
case of the density operator) and its Fourier transform (FT). Thus, a unitary operator U 
corresponding to a given surface in double phase space is represented by U(Q) = {q+\U\q-), 
U{P) = {p+\U\pJ), or in the Weyl representation as U{q,p), depending on rotations in double 
phase space that correspond to different choices of FT’s in quantum mechanics. 

But the analogy of an operator represented by a lagrangian surface in double phase space 
to the state represented by a lagrangian surface in simple phase space can now be pushed a 
step further: What about the ’double Weyl transform’ of a superoperator, represented as a 
function in double phase space? The answer invokes the adaptation of the Choi-Jamiolkowsky 
isomorphism [|5l, well known in the theory of quantum information, to the operator basis 
underlying the Weyl representation and its Fourier transform. The purpose of this paper is 
to undertake this adaptation, showing that the Choi or dynamical matrix of a superoperator 
is the ’double Weyl transform’ of its matrix elements in the Weyl basis. We achieve this by 
defining reflection and translation superoperators, labeled by points in double phase space, in 
strict analogy to the corresponding definitions in single phase space. Moreover we show that 
they take a simple monomial form in the Choi conjugate basis. 

Another surprise is that the Choi-conjugate pairs of operators representing pure states 
are identical. This leads to unexpected Fourier identities involving integrals of products of 
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ordinary Wigner functions. In some cases this in turn leads to new identities involving the 
special functions of analysis. 

The paper is organized as follows: in section 2 we review the concept of double phase 
space [HI El El [H and how it is related to the labeling of conjugate bases of operators in terms 
of lagrangian surfaces in double phase space. 

In section 3 we extend the methods to the representation of superoperators. The concept 
of Choi-conjugate bases is defined and applied to the phase space bases that are determined 
by reflections and translations. In this context Choi-conjugation takes on the character 
of a “double” Wigner -Weyl transformation in double phase space. The identification of 
such double Wigner functions with the Wigner-Weyl transformation from alternative matrix 
representations of a superoperator is examined in section 4. 

Superoperators for various kinds of quantum evolution are then discussed in section 5. 
Finally, in section 6 the pullback of the double phase space results for single phase space 
Wigner functions generates several new identities for pure states. 

2. Double phase space and operator representations 

The concept of double phase space underlies in classical mechanics the general theory of 
generating functions of canonical transformations. It provides an elegant way to visualize a 
canonical transformation in a doubled phase space as the gradient of a generating function 
HI El El- In quantum mechanics it provides a flexible mechanism to represent unitary 
propagators in the semiclassical limit IfTOll as functions on these manifolds and underlies the 
presentation of quantum mechanics in terms of phase space path integrals liTTlfT^ . 

A canonical transformation [f] x_ i-^ x+ can be specified implicitly by a generating 
function S {q+, q-) whose differential is 

d5 (<?+, <?_) = p+d<2'+ - p_dq_ (1) 

With a simple redefinition of coordinates 

2i = {q+,q-) Pi = ip+,-p-) (2) 

we rewrite d5 = Pi ■ AQi and the transformation is defined implicitly as PfQi) = 
dS{Qi)ldQi. The new coordinates Qi,Pi can now be interpreted as canonical coordinates 
in a phase space with doubled dimensions. In the elementary theory other canonical pairs 
are well known and are obtained by various Legendre transforms. In the standard notation of 
1131 HU [TSl they are: 

ds = dFi = p+dq+ - p-dq. = Pi ■ dQi Qi = (<?+, Pi = {p+, -p_) 

dF2 = q+dp+ + p.dq. = P2 ■ dQ2 Q2 = ip+, <?-), Pi = iq+,P-) 

dFs = - q-dp- - p+dq+ = -Q 2 - dP 2 

dF 4 = - q.dp. + q+dp+ = -Qi ■ dPi (3) 


f For simplicity from now on we restrict to the 1-D case x = (q, p) 
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The general theory developed by [[B IS 13 interprets these alternative pairs Qj, Pj as the 
eoordinates for hyperplanes lying in a phase spaee with doubled dimension. We give here a 
synthetie overview with the purpose of identifying the structures that will reappear in quantum 
mechanics in the representation of operators and superoperators. Starting from the elementary 
statement of area preservation in a canonical transformation 




p.dq. 


rewritten as 


- q; x+ ■ 7dx+ - x_ • 7dx_ = 0, 


(4) 

(5) 


where as usual x+ = (q, p)+ are points in single phase space and ^ 0 j standard 

symplectic form in 2D. Placing ourselves in the 4D direct product space spanned by (x+,x_) 
we rewrite Q as 



( 6 ) 


Thus the canonical transformation x_ x+ is defined by a 2D lagrangian surface in the 4D 
direct product phase space. On this surface 



(7) 


is an exact differential. However this surface is lagrangian with respect to the non standard 
symplectic form | ^ j. To introduce canonical coordinates with respect to the standard 

canonical form ( ^ |, we consider a linear transformation 



such that 

dA = ^(x+ ■ 7dx+ - x_ ■ 7dx_) = dQ-Q- dP). 
The matrix Id should satisfy 




J 

0 



0 

1 


-1 

0 


( 8 ) 

(9) 

( 10 ) 


Various choices of the matrix 'I/, lead to different pairs Qt, Pi, and all such pairs are related 
by canonical transformations in 4-D phase space. Notice that the first form of dA in Q 
is invariant under canonical transformations x+ ^ y+ in single phase space while the 
second form becomes invariant under full two degrees of freedom canonical transformations 
Q,P Q',P' ■ To proceed to the definition of the canonical transformation in these 
coordinates a Legendre transformation on dA leads to 

dA = d(A + ^P-Q) = P- dQ 


( 11 ) 
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P{Q) = 


dA{Q) 
dQ ■ 


The eanonical transformation is then speeified via the parametrie equations 


/ x^(0 \ / t/oo f/oi \/ Q \ 

\x-( 0 / I 00 01/100/ 


( 12 ) 


(13) 


Alternatively a different Legendre transform leads to a eonjugate representation exchanging 
coordinates and momenta dA(/’) = d(A - }^P ■ Q) = -Q ■ dP leading io Q = -dAjdP and to 
the parametric equations 


/x,(0\ /t/oo 0l\/00\ 

\ X_(0 / \ 00 01 /\ Z’ / 


(14) 


In this general approach a canonical transformation is characterized by a) a linear 
transformation 14 satisfying (10) and defining the canonical coordinates Q,P and b) by 
generating functions A{Q) iA{P)) on the lagrangian surfaces P = const. {Q = const.). The 
two generating functions are themselves related by the Legendre transform A = A + P ■ Q. 
Singularities in the process of inverting the functions x_(2) (x+(/’)) in ( [T3] ) and ( [T4| ) will lead 
to transversality conditions guaranteeing the existence of the transformation. This general 
scheme can be easily identified in (|^ where Q\,Pi and 0,0 are canonically conjugate pairs 
and '1/ is a simple permutation. 

A different choice of canonical variables in double phase space is related to the well 
known Wigner-Weyl representation of quantum mechanics in phase space. It is characterized 
by the matrix 




-U 


leading to the transformation 

x+ = 2 + ^JP 

and its inverse 


x_ 


Q--,JP 


(15) 


(16) 


2=^^^^ P = -Jix,-x_). (17) 

The lagrangian planes of constant Q represent the transformation x+ = 22 - x_ which is a 
phase space reflection, while those of constant P lead to the translation x+ = x_ + JP. Thus 
Q,P are natural labels in quantum mechanics for the operators that represent translations 
and reflections. However it is customary to use, instead of (Q,P), the single phase space 
coordinates 


X = 2 = (x++ x_)/2 ^ = JP = x+-x (18) 

They have a simple geometrical interpretation in single phase space as centers and chords of a 
pair of phase space points ifT^ and the Poincare generating function, S (x) has many desirable 
properties [|9l[T^. 
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Lagrangian coordinates in double phase spaee provide the natural labels of operator 
bases in quantum meehanies. In the simple ease of the position basis \q) an operator is 
represented by its matrix elements {qf\d\q-) = ix d\q-){q+\ as a funetion on the lagrangian 
plane Qi = (^+, q-). In the conjugate momentum basis \pf){-p-\ the same is true as a funetion 
on Pi = ip+,p-). 

To treat the general ease we eonsider a basis of operators S „ labeled by a double index 
or 6 whieh ean be any of the lagrangian planes eonsidered. Introdueing the notation 
((A|ifi)) = txA^B for the Hilbert Schmidt scalar product, we adopt a dual notation for operators: 
when they are considered as veetors in the Hilbert-Schmidt sense we use the double Dirac 
notation ||A)), but when we eonsider them as two sided arrays with the standard multiplieation 
rules they will be denoted simply as A. The latter notation emphasizes the aetive aspeet of 
A as operators, while the former is meant to bring out their passive role as basis elements. 
The bases eorresponding to the lagrangian planes in Q are \\Qi)) = \q+){q-\, ||Pi)) = 
|p+)(-p_l, 1122 )) = \p+){Q-1 11 ^ 2 )) = \q+){P-\- Their eomplementarity is refleeted in the 
property ((P,j|2,)) = exp -iQi ■ Pi/h. The eommon properties of these bases are 


«5J|54 = Ad(^r-;8) 
i J d2a||S„MSJi = l 


Orthogonality 

Completeness 


(19) 

( 20 ) 


One should note that, within the context of the following seetion, the latter relation aetually 
represents the identity superoperator. The faetor A is eonventional and could be removed by 


a simple resealing. Eq.(20) implies that any operator O ean be expanded in the basis ||5 q.)). 


i.e. 


0=1 f CISC'S AS JO}}, (21) 

with eoeflieients ((5a||0)) that are the e-number representatives of O . In the simple ket, bra 
bases eonsidered above they are simply its matrix elements. 

Another pair of complementary operator bases is obtained when considering the unitary 
representation of refleetions and translations in phase spaee. They are defined in terms of 
symmetrized Fourier transforms of the previous bases as follows 

/ / t 

dc{\q^%^}{q-'^-\A‘i ( 22 ) 

^ = rdp' Ip' + |)(p' - || (23) 


Here the first is in the position basis and the latter in the momentum basis, but either ean 
be used. [§] To simplify notation we have defined to = exp i//i. As a complementary set of 
operator bases they have the properties 


{{RJRA = 2;r^d(x - x'), {{TJTA} = 27rM(^ - f), ITJRJ) = ^>“<^’*>(24) 


so that for these bases A = 2nh, whereas A = 1 for the position or momentum bases. The 
action of these operators and further properties are reviewed in the Appendix. 


§ It should be noted that the unitary reflection operator is Px/2, but our deflnition is more convenient as a basis. 
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The c-number representatives of operators in these bases are phase space functions 
that constitute the starting point for the formulation of quantum mechanics in phase space 
[fT71 [m [13 [20l . The reflection basis, 


6>(x) = {{Rf\d)) = 



+ ^\0\q- |-) m 


■q'p 


(25) 


leads to the Wigner- Weyl symbol of the operator while the translation basis 


0{^) ^ iT^WO)) = Jdp' ip' + ||0|p' - |) (26) 

leads to its characteristic function. In line with [fT^ we will refer to them respectively as 
centre and chord representations. They are related by a symplectic Fourier transform 


0(x) 


-2^nf 


■2 n 


dr 


Oi€). 


(27) 


In the special case of the density operator p, p(x) and p{f) are the Wigner function and the 
characterisitc function, or chord function: 


T(^) 


1 


2nh 


«rilp)) and 


with the special normalization properties 
InhxiO) = trp = 1 and 
corresponding to classical distributions. 


W(x) 


1 


2nh 


m\p)). 


f 


dx^ VF(x) = trp = 1, 


(28) 


(29) 


3. Superoperator representations 

3.1. General representations and the Choi- Jamiolkowsky isomorphism 

A superoperator S is a general linear transformation in the space of operators. Simple 
examples are the unitary evolution of a density matrix p ^ UpU\ symmetry operations, 
time reversal, non-unitary evolution, etc. Once the action of S on general operators has been 
specified as S(0) = O', its matrix elements in the S a basis are 

{{SJS\\Sp)) = tx[SlS{Sfi)] (30) 

Thus a general superoperator can be expanded as 

S = ^ J d^ad^fi «S Ji S ||5^))||5„))«5^|i. (31) 

The superoperators ||5 'q,))((5j 8|| constitute then a complete orthonormal basis labeled by a 
double index a,P e Their action on an operator O follows the standard bra and ket 
rules so that \\Sa)){{S^\\0)) = S^trS^O. The unitary relations between different choices of 
bases also allow us to define the trace of a superoperator uniquely as 

TrS^^Jd^ct «5 Ji S ||5„)). 


(32) 
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Nonetheless, it is eonvenient to eonsider a different aetion related to the pair of operators 


j. 

c cT 


{0) = SaOS^ 


(33) 


Indeed, S is a different orthonormal superoperator basis related to the former by a unitary 
transformation in The transformation is given by 

= f tr(Sj^S;sl) Ii5,))«5^||, (34) 

with the inverse: 


l|5a))((5^|| 




(35) 


We prove for example (34) by eomputing with (30) the matrix element ((5 q.|| • 5], HSyg)) = 

tr(SxS pSyS 1). In the following we will refer to the pair of superoperators ||A))((5|| and A* 
as Choi-eonjugates and extend this name to the two bases. 

S ean also be expanded in the new basis as 


S = ^Jd^xd^yCs(x,y)S,, • S^, (36) 

so that the coefficients Cs(x, y) constitute an alternative representation of S, related to its 
matrix elements by ( [35] ) 

«5 J| S 1154 = ^f dWyCs(x,y) triSlsJ^S^). (37) 

The nature of the transformation between the two representations of S is specific of a 
given basis and depends on the trace of four basis elements. The simplest transformation 
occurs when a transition basis E/j = |/)( 7 | of orthonormal states is used. In that case (setting 
A = 1) the quadruple trace results in a product of four delta functions which reduce the 
relationship to a partial transposition of indices 


((E,jHSmd} = Cs(ik,j7). (38) 

In this case one obtains two different linear transformations related to S. They are well known 
in the quantum open systems literature. Csiik, jl) is known as the Choi matrix ETl . dynamical 
matrix [|2^|23l, chi matrix [|24l and the relationship between them as the Choi- Jamiolkowski 
isomorphism. It is also of importance in the literature on matrix product states for spin chains 
ll25ll . When S represents a quantum operation on density matrices the physical requirement 
of complete positivity requires that Cs(ik, jl) be hermitian and positive [[2^ |2T]|. In that case 
Cs maps the space of positive density matrices into itself. Although related, the spectral 
properties of the two transformations are very different. Thus if Cs(ik, Jl) is hermitian it can 
be brought to diagonal form leading to the Kraus representation Il2^ . On the other hand 
{{EijW S \\Ei,i)) as a linear transformation may not even be normal, leading in general to a 
complex spectral decomposition. 
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Consider first the refleetion basis ^x- A general superoperator S in this basis has the matrix 
elements S ||^x_))- The labels, x+ in R^, span a full set of refleetion eentres, eaeh 

covering a full phase space. On the other hand a rotation in double phase space allows us 
to consider ||i?x+))((^x_ll as an alternative position basis for superoperators, just as valid as 
\qf){q-\ in the operator case. Hence, one defines a double Weyl transform in analogy to ( [25] ) 
as a symmetrized Fourier transform [jl] 




<Xl,^> 


(39) 


where we use the phase space coordinates (x,^) which change the Fourier phase from -Q ■ P 
to Xi • =< Xj,^ >. This leads to a very natural definition of a reflection superoperator in 

this basis as 


Rx,^ = 


— r 

2nhj ‘ 


-<X|,^> 


d Xl I|i?x4x,))((^x-ix,ll 

The Weyl transform of S can now be expressed as a superoperator trace, 
S(x,^) = TrSRx,f, 


(40) 


(41) 


in complete analogy to (25). The “super” reflection properties of Rx,^ are contained in the 
easily derived actions on translations and reflections, 

Rx,^ll^xo)) = 


Rx,flirfo)) = 


(42) 


in obvious similarity to (A.6). 

Just as any other superoperator, one can now expand Rx^ in the Choi conjugate basis, 
Rxi • Rx 2 , using the general relationship in (35): 


||J?x.ix.))((^x-ixJI = 


{Inh) 




Ra • Rp ir{RaR 


X+ixi 


RpR 


lx,)- (43) 


X-^Xi 


The quadruple trace computed in (A.13|) reduces this to 


||J?x.ix,))((^x-Ax,ll 


a -V P 


= ^ Jd^rz d^P K • Rp d(x - 




and its inverse 


^x+f/2 • Rx-^/2 


- f 

2nh J 




<Xl,f> 


d^xi ijR^^,^f)((R^iJioj- 


-<X|,f> 


(44) 


(45) 


Comparison with (40) then shows that Rx,f has a very simple (monomial) form in the Choi- 
conjugate basis: 


Rx,f - R 


x+f/2 • Rx-^/2 


= Rx •Rx 


( 46 ) 


II The denominator 2nfi = A is usually absent because A = 1 for the position representation. 
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(47) 


Very similar results ean be obtained in the translation basis. A translation superoperator 
is defined as 

^ ^ a.-'*--, 

with the translation properties 

Tx,fll^xo)) = ll^xo+x)) 

It defines the double ehord representation of S as 

An almost identical derivation as in ([43]) yields the relationship between the Choi-conjugate 


TrSX:^,. 


(48) 


(49) 


translation bases: 

A A + 

T.+lr • 

*^2® X-if 




<^1,X> 


=- f 

Inh J 


which again lead to a monomial representation of the translation superoperator: 








A A i 

Tx 

x+ ^ X- 


(50) 


(51) 


Equations (46),(51) are our main general results: in double phase space the reflection and 
translation superoperators can be defined as usual in terms of symmetrized Fourier transforms 
of matrix elements, but they also aquire a simple monomial form in their respective Choi- 
conjugate basis. The interplay between their action in double phase space with the single 
phase space factor translations and reflections in Rx^»Rx and rx^.»rx_ will be further discussed 
in section 5. 

We now use these identities to derive the Choi-Jamiolkowsky relationship between this 
pair of conjugate bases that represent S. Expanding a general superoperator as in ( j^ 

d^x+d^x 




{2nny 


Cs(x+,x_)i?x^ •Rx 




d^x+d^x. 

(2nhf 


Cs(x+,x_)rx^ • fl. 


(52) 


defines Cs(x+, x_) and Cs(x+, x_) as the Choi or dynamical matrices of S in the reflection and 
translation basis. Changing the double phase space variables from (x+, x_) ^ (x,^) as in ( |T8] ) 
(with unit jacobian) we rewrite it as 


-f 


d^xd^^ ^ f 

^Cs(x+-,x--)R.^ 


-f 


d^xd^^ ~ ^ . 

-|Cs(x-h |,x- -)T 


x,^, 


(53) 


(27rhf 2 2' 

that is, the superoperator is expanded as a superposition of reflections (translations), such that 
the coefficients are identified with the center (chord) representations in double phase space. 
Thus, we have obtained 


S(x,^) = Tr(SRx,^) = Cs(x^,x_) 
S(x,^) = Tr(STl^) = Cs(x^,x_), 


(54) 

(55) 















Representation of superoperators in double phase space 


11 


which is the desired Choi-Jamiolkowsky relation in these bases: the Choi matrix appears as a 
rotated double Weyl or ehord transform of the superoperator. The property ean be inverted to 
retrieve the matrix elements from an integral over the Choi matrix 


{{R 


X+ixi 


S l|i?x-lx.)) 




=- r 

Inh J 

= - f 

Inh J 


d^^ Cs(x+ -^,x- -Ooj- 
d^x Cs(x+ ^^,x-^^)co 


-<X|,f> 


-<^|,X> 


(56) 


It is important to note that these relationships follow direetly from (37), but the direet 
derivation would then miss the interpretation of the Choi matrix as a phase spaee distribution 
in double phase spaee. Furthermore, the complete analogy between the Weyl transform in 
double phase spaee to its familiar version in ordinary phase spaee allows us to direetly import 
formulae that are familiar in the latter eontext into the superoperator seenario. For instanee, 
the traee of a superoperator, S, takes on the alternative forms: 


Tr S 


=- r 

iTTh J 

= - f 

Inh J 


d^x m\ s \\R,)) 


{2nh) 


kJ 


dhd^^S(x,i) 


d2^«r^lisiir^)) = s(x = 0,^ = 0). 


(57) 


Again, the traee of a produet of superoperators may be expressed as 


Tr S2 Si 


/ 

/ 


d^x.d^x^ 

{2nhY 


((^xJI S2 ||i?x,» ((^xjl Si \\R,_)) 


I 


|2„j2v 


(2^^)2 HT^Jl S2 ((r^JI Si ||r^_)) = J dW^S2(x,^)Si (x,^). (58) 

The Wigner-Weyl expressions of partial traces over subsystems ean also be immediately 
generalized from their versions in single phase spaee, whieh are discussed in [[81. 


4. Generality of the double Weyl transform 


In the previous seetion the double Weyl transform was implemented by ehoosing the eonjugate 
lagrangian planes Q, P (or x,^) in ( [T^ . We now show that many other equivalent definitions 
ean be defined using alternative eoordinates. Take for example the eoordinates Qi,Pi in ([^ 
and the operator basis 

Qu = WQJ = Pa = \\Pa)) = \p^){-P-\, (59) 


a = {q+,qY,0! = (p+, -P-) are eonjugate position and momentum variables in double phase 
spaee. Starting again from the double position basis ((2allS||2a')) an alternative double Weyl 
transform in this basis is obtained in exaet analogy to ( [39| ) as 

S(a,a) = JdV << 23 , 13 ,II S||<2a-.a'))^^“"' (60) 

Again a refleetion superoperator results as 


R'(a,or) 


/ 


dV || 23 .ia'))((ea- 


a'l 


to 


(61) 
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where now the refleetion is in the position basis \\Qa))- In this ease the Choi eonjugation ( [35] ) 
between the two position bases eonsists of a simple transposition of indiees 


\\Qa,b)){{Qc,d\\ - Qa,c • Qb,d 

Applying this transformation to ( |M] ) 


(62) 


^+11 d +, 


R'(a,^r) = J dq',dq'_\q, + • \q_ - (63) 

we reeognize again the appearenee of the refleetion superoperator -Thus 

R'(a,a) = ^x,«^x_ (64) 

where (a, or) = {q+,q-, p+,-p-) is related to (x+,x_) = {q+,p+,q-, p^) by a signed 
permutation matrix ddi that satisfies ( [T0| ). Just as in (54) we then obtain again 

S'(a,or) = Cs(x+,x_). (65) 


Thus a rotation of the arguments of the Choi matrix (by a linear transformation 'l/i satisfying 
( [T^ ) results in the double Wigner transform of S in a different basis. Of eourse, the reversal 
of the Weyl transform ( |5^ ean also be used to retrieve the double position or the double 
momentum representation: 


HQa 


a'li S IIQa-ia')) = f S'(a,«) 


Moreover, eomparing ([65]) with (|54|) we also obtain 


( 66 ) 


S'(a,or) = S(x,^) = Cs(x+,x_). 


(67) 


In short, we obtain that the single function Cs(x+,x_) supports the Weyl transform of S in 
two different bases just as would be expected from symplectic invariance for the ordinary 
Weyl transform in single phase space. We note here that the classical transformation (a, a) 

(x, -J^) preserves the symplectic form in 4D and therefore belongs to the symplectic group 
Sp(4). An almost identical derivation involving the bases II22)) = \p+){q-\,\\P 2 )) = \q+){p-\ 
yields S"(Q 2 ,P 2 ) = Cs(x+,x_). Clearly an analogous discussion relates the various double 
chord transforms to the Choi matrix Cs(x+,x_) in the translation basis. 

Indeed, the immediate generalization of the above discussion would show that the 
transit between the various operator representations corresponding to different choices of 
double phase space coordinates (|^ among themselves (as well as with the Weyl and the 
chord representations) is universally obtained from a single double Weyl representation by a 
symplectic change of coordinates. To complete the generalization one would need to attempt 
the explicit construction of metaplectic superoperators in double phase space that would 
implement the change of basis. We intend to proceed along this path in future publications. 


5. Quantum evolutions 


Consider the superoperator 

v = u*u\ 


( 68 ) 
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which propagates unitarily the density matrix as Up Uf The matrix elements in the refleetion 
basis are U ||i?x_))- This is the kernel that propagates unitarily the Weyl representation 
of the density matrix {{RxWp)) 






d2x_ {{RxJ\ U \\RxMRxJ\p)). 


(69) 


As a particular case of the operators treated in section 3, the Choi matrix of U is easily 
eomputed in a separable form as Cu(x+,x_) = t/(x+)t/*(x_) where U(x) = {{Rx\\U)) is the 
Weyl transform of U, i.e. the Weyl propagator. The double Weyl transform of U is then given 
by (60) as U(x,^) = TrURx,^ = U{x+ ^^)U*(x - The double propagator ean then be 
computed explieitly using (|5^: 


U ^ Jd'f U(x + (70) 

At first sight, this eentre-eentre propagator seems to result from a double Weyl transform from 
the single Weyl propagator, U(x), of the unitary map U (see [12711^), but the derivation of 
( [5^ clarifies its true role as the inverse transform from the double Wigner funetion, that is the 
Choi representation of the super evolution operator. ^ In the semielassieal regime the ordinary 
Weyl propagators have explieit formulae in terms of generating funetions OTl IT^ . Further 
evaluation by stationary phase leads to semielassieal approximations for this propagator. Very 
similar formulae for the chord-chord propagator ean be derived in the translation basis: 


((T^JI u ||r^_)). 

A more general evolution is generated by the Kraus superoperator 


K = (71) 

j 

which is again easily aeeommodated in the present framework. Indeed, linearity of the Fourier 
transforms then specifies the integral kernel for evolving Wigner functions as 

«/;„.„||K|K_,„» = J d^#/f,(x + if)Jf/(x-if)£j--f>. (72) 

An alternative generalization from ordinary unitary evolution, whieh is ineluded in the 
general formulae of the previous section, is the superoperator 


U = Cl • ul- 


(73) 


This determines the evolving kernel for the quantum fidelity, or the quantum Losehmidt eeho 
||^I^I3?1 , that is, the overlap of two different evolutions for the same initial state. Unlike 
the previous examples, the trace of the evolving operator is not preserved and the identity 
operator is not invariant. 

The super-refieetion, Rx,^, in (46) is preeisely of this form (within faetors of two), 
ehoosing Ui = R^+y and U 2 = Rx-\^- The same goes for the super-translation, Tx,|:, defined in 
( |5T] ). As revealed by ( [42] ) and ( [48] ), these superoperators can be viewed as active agents, rather 
than mere passive Choi bases. For instance, the general rules for produets of refleetions and 


K One should note that this product of simple Weyl propagators in the integrand was mistakenly identified in 
Ii30l with the mixed centre-chord propagator, defined in Ii29l . 
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Figure 1. The symplectic area of the quadrilateral, A 4 (^i,^ 2 ’^ 3 >^ 4 ) - 2 (< > 

+ < ^ 3,^4 >), determines the phase for the action of the superoperator and the 

superoperator ^^4 • ^X 2 ■ It is also identihed with the phases for the trace of the product of four 
reflections or four translations that are evaluated in the Appendix. 


translations, which are reviewed in the Appendix, determine the aetion of a super-translation, 
aeting on the translation operator as just with^j = ^ 2 +^ 3+^4 as depieted 

in Figure Likewise, the aetion of the superoperator ^X 4 • ^X 2 oa the refleetion is simply 
with Xi = X 2 = X 4 - X 3 . Note that the four reflection centres Xy form a parallelogram 
with half the sympleetie area of the cireumseribed quadrilateral, A 4 , formed by the translation 
ehords, ^ j, so one obtains the same phase for the eomposition of three translations or three 
refleetions. Indeed, this is also the same phase as determines the traees of four translations or 
four refleetions that are diseussed in the Appendix. 

The extent to whieh ( |4^ and ( [5T] ) are natural definitions for the reflection and translation 
superoperators is iluminated by considering the graphs of classieal translations and reflections 
as planes in double phase spaee. First one notes that eaeh operator, eorresponds to the 
translation, x_ x+ = x_ -1-^3: The loeus of all pairs (x_, x+) in a given initial plane in double 
phase spaee. Then the pair of translations, x+ 1 -^ x^ = x+ -l- ^4 and x_ x^ = x_ - ^ 2 , that is, 
the double phase spaee translation X X', takes all double phase spaee points in the original 
plane to the new plane x+ = x_ . Thus, the new and the old planes representing translations 
are themselves related by just the double phase spaee translation eorresponding 
that is, = to In the same way, one also finds that the eombination of phase 

spaee refleetions within eaeh of the phase spaces x+, eorresponding to R^^ • R^^, transport the 
double phase spaee points X = (x_ = X 3 - ^^ 3 , x+ = X 3 - 1 - ^^ 3 ) in an initial refleetion, to a new 
reflection, Xi -l- 1 -^ Xj - 


6. Phase space correlations and general pure state conditions 

Consider a general quantum state represented by its density matrix p. Its phase spaee 
translation eoresponds to the operator 

A. A J. 

A rri A T 

Px = T^pTf 

and the (non-normalised) phase space correlation function lIMl fSlI is naturaly defined as 


( 74 ) 
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Likewise, from the refleeted state = R^p Rx, an analogous anticorrelation function can be 
defined as 

C’‘ = trpp’‘ = tr(pRxpRx). (76) 

To provide a link with our previous results it is then eonvenient to define the density 
superoperator, 

^=P*P, (77) 

in terms of whieh one can identify 

Cx = iTx\\^\\Tx)) and = {{Rx\\V \\Rf)) (78) 


as the diagonal matrix elements of P in the translation and reflection bases. Utilizing the 
cyclic invariance of the trace in ( [75] ) and (761, we can alternatively rewrite 

Cx = {{p\\Tx*Tl\\p)) and = {{p\\Rx* RxWp))- (79) 


The Choi-conjugate relations ( |50| ) ( |45| ) lead to 


jXl 


lip)), 


(80) 


( 81 ) 


that is, both correlations can then be computed directly from the chord and the (ordinary) 
Wigner function and chord function ( [28] ). 

In the case of pure states, p = |<A)(iAl, the correlations can be computed directly: 

Cx = |((rx|ip))|2 = {2nhf\x{xt and C" = |((i?xllp))P = {2nhfW{xf (82) 

and then comparison with ([80l),([M]) leads to the identities [|35ll36l, 

1 t(x)I^ 


- r 

2nh J 




and 


1 T(x )2 


/ 


d^v y y 

^lT(x.|)lT(x-|). 


(83) 


(84) 


Further relationships between products of center and chord functions for pure states are 
obtained from the off-diagonal matrix elements of P, such as the matrix element 

«r,-ifll P = ((pll J-„i{ • Tl iw 


^ f (85) 


2nh 

Now we introduce the remarkable property of pure states that the superoperator P is self 
Choi-conjugate, i.e. p • p = ||p))((p||, as can be easily checked by acting on any operator. 
It is then possible to pull back the Choi-conjugation relations as unsuspected properties of 
ordinary pure state Wigner or chord functions in single phase space. Indeed the l.h.s. of (851 
is then reduced to ((r,j_i^||p))«p||r,j^i^)) so that 


;r(x + - W = ^f T(^i + -20r(^i - 


( 86 ) 















Representation of superoperators in double phase space 


16 


is identified as a general requirement for chord functions of pure states. The previous case 
of (831 is included as ^ = 0. The off-diagonal elements in the reflection basis yield a similar 
property for Wigner functions, 

W(x + W(x = 

and a mixed case relating Wigner and chord functions is 

W(x + i^) W(x - i^) = ^ Jd^^i 

These formulae establish families of Fourier identities parametrized continuously by 
X or They show e.g. that the product of symmetricaly displaced Wigner functions, 
VF(x -I- 1^) W(x - is its own Fourier transform for all x and Xi^ + is 

likewise invariant for all values of Many particular cases give rise to interesting properties. 
Besides the ones already noted for the correlations, we list the following special cases 


(87) 


( 88 ) 


XiW = ^f 

W(i^) W{-\^) = ^f d\ W(ixi) W(-7xi) m 




(89) 

(90) 

(91) 

(92) 


W(i^) W(-i^) = ^ J X(^ + yi)x(^ - -2^i) 

i^(x)" = ^ J d^^:X(^^i)X*H^i) 

The last one, taking into account that;\f(x) = ;^:'*(-x) and Parseval’s relation, yields the further 
integral 

Jdh W(x)" = Jdh [Y(ix)r. (93) 

Furthermore, for x = ^ = 0, one obtains: 

W^(O)' = ^ Jd^x W(ix)W(-ix) = ^f d^^,x(^^i)xH^i)- (94) 

The way the general formulae ([8^,([87]),([88]) were derived implies that they are necessary 
conditions for pure state distributions. By setting x = 2^, it is easily shown that (86) is 
equivalent to the pure state condition p = in the chord representation, and so it is also 
sufficient. The case of the Wigner function is not so transparent, but the integral product rule 
for the Weyl representation [flbl for p^. 


p2(x) = 4 
= 4 


f 

is imediately simplified by ( [87] ) so that 

J 


d^Xid^X2 W(xi) W(X2) CO 


2<(x|-x),(x2-x)> 


d^^ W(x+i^) W(x-i^) m 


2<(x-x),f> 


p (x) = STth VF(x) d X VF(2x - x) = IttH VF(x) = p(x). 


(95) 


(96) 
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There are some notable eases where the Wigner funetion of pure states ean be deseribed in 
terms of standard speeial funetions found in eg OTl . Then the identity, 


lT(x+|) lT(x-|) 


f 




(97) 


implies a possibly unsuspeeted Fourier identity for a symmetrized produet of analytic 
functions. Such is the case for the eigenfunctions of the harmonic oscillator - given in terms of 
Laguerre polynomials IfT^ - or the unnormalized eigenfunctions of the hyperbolic hamiltonian 
H(x) = p q, calculated in OSlI in terms of Laguerre functions of complex index. We develop 
here the important example of the linear potential V{q) = q. If m = 1/2 and h = 1, so that 
the Hamiltonian is simply H(x) = p^ + q, the zero energy eigenfunction is proportional to the 
Airy function, 

r / r,3 

which is not normalizable and hence has no Fourier transform. Evidently, the corresponding 
momentum representation of this state is just 


dp exp 


1 y + /^<? 


(98) 


(p|(A) = 


1 


^/2n 


exp 


. P 


(99) 


so that the corresponding Wigner function is just lEUI 


W(x) 




(p + - ^)exp[//^] 


exp 


,'3 


i\^ + ip^ + q)p' 


21/3 


n 


Ai [l^l^H{x)), 


( 100 ) 


which is also not square-integrable. 

Nonetheless, the product VF(x + y) VF(x - y) decays exponentially in the x'-phase plane 
outside the region limited by the pair of reflected parabolae, H{x + y) = 0, so that it is square 
integrable. Thus the direct verification of the Fourier identity ( [87| ) in the case of the Airy 
function proceeds from the integral representation: 

(2^/rfp'rfp-'exp 


^ + (((>1 +f)"+9i + f IP 


exp 


.,//3 


.,P L P2,2^ d2. , 

-M —-|(/^i-y) +^i-yb 


12 


( 101 ) 


Then the transformation p' = a + p" = a 


I simplifies this into 


VFixi + ^iWlxi 


2 


( 2 ^/ 


da db exp 


b'^ a^b 


P2 


48 


-I—-—h {q 2 + 2pip2)a + (pi -\ —-—I- qi)b 


.( 102 ) 
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Figure 2. The product of two Airy functions as in 
density plot blue is positive and yellow negative 


( 104| ifor X 3 for Xi 


= (-3,0). In the color 


Thus, the Fourier transform with respeet to X 2 beeomes 


^ J'dx2 ^+ y) W" (xi - y) = J'da db dq2dp2 


exp 


., b^ a^b N / 2 P'i 

1 1 ^ ^ + (<?2 + 2piP2)a + (pi + — + qi)b + p2<?3 - diPs 




db dp 2 exp 


b^ , Pi^b , ^_ , ^ 2 , P2 


1 I - + 

48 4 


+ 2piP2P3 + (Pi +^+qi)b + P2q3 




(103) 


which implies the Fourier invariance for a product of symmetrized Airy functions, 

^ Jdx2 + ^)) Ai(22/3//(xi - ^)) 

= Ai ( 2 ^ 1 ^Hix, + ^)) Ai ( 2 ^ 1 ^Hix, - ^)), (104) 


not encountered even in [|39l , a book dedicated specifically to Airy functions or [|40l . In 
Figurewe show this product in the plane X 3 for Xi = (-3,0). Notice the central symmetry 
implied by the Fourier invariance. 

An important feature of this example is that the semiclassical transitional approximation 
of pure state Wigner functions for general WKB-quantized states near the closed energy 
eigencurve was shown by Berry [|4T]| to be just the Airy function over the approximating 
parabola. Thus one finds that both the exact Wigner function and its transitional 
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approximation satisfy the new Fourier invarianee, even though this is not the ease of other 
semielassieal approximations that are more refined in other respeets. 

7. Conclusions and outlook 

Two possible representations of a superoperator can be naturally derived from the same 
operator basis. They are unitarily related, and we have referred to them as Choi-conjugate. 
We have developed here the general relationship between them - when the basis is orthogonal 
- and we have studied in particular the form of this relationship when the unitary bases of 
translations and reflections are used. It turns out that the representation in terms of the Choi 
or dynamical matrix Cs can be interpreted as a double Weyl or Wigner transform of the matrix 
elements of the superoperator. This is because the representation in terms of the Choi matrix 
is actually an expansion in terms of translation and reflection superoperators, in strict analogy 
to the expansion of an ordinary operator in terms of translations and reflections, yielding its 
Weyl or Wigner transform. The definition of Tx,|: and opens up the possibility for a full 
study of the affine geometry of these superoperators in double phase space, including the 
definition of unitary superoperators that implement symplectic transformations belonging to 
S p(4). We intend to pursue this analysis in the future. 

Our treatment here has been for the simplest case of a phase space with no boundaries and 
with one degree of freedom. The extension to D degrees is immediate and needs no further 
comment. The adaptation of our techniques to a phase space with boundaries needs more care. 
The case of torus topology - periodic boundary conditions both in position and momentum 
-is the closest to the present approach and leads to a finite dimensional Hilbert space of 
integer dimension d = area./2nh [I^ . fl4^ . [l44l . ll45l and is of great current interest in quantum 
information theory. Translation and reflection operators can still be defined and provide a 
basis for a similar treatment as the one developed here. In this context, the display in double 
phase space of the properties of superoperators can provide new insights into their actions, 
just as the celebrated Wigner and Weyl representation displayed properties of quantum states 
in single phase space. The action of gaussian noise channels in the chord representation Il4^ 
is a first step in that direction. 

As application of these methods we have found some previously unknown identities 
relating products of Wigner and Weyl distributions for pure states. These identities generalize 
the pure state conditions and in some cases produce new relationships for the special functions 
of analysis. 
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We review the well known definitions and properties of refleetion and translation operators. 
They eonstitute the foundation for the Weyl representation of quantum meehanieal operators 
as phase spaee e-number funetions (the Wigner quasiprobability distribution in the ease 
density matriees). We start with the usual q,p operators that we subsume in a phase space 
operator x = {q, p) and a phase space label x = {q, p) 6 The corresponding position and 
momentum bases are denoted \q), \p). The symplectic product is defined as 


< x,x' >= {q,p) 


0 -1 

1 0 


q 

P' 


pq -qp . 


For notational simplicity we also introduce the quantities 


T = 


m = e'^l 


(A.l) 

(A.2) 


Reflections and translation operators are defined as 

/ t f 

\q + \'i^q- yl ^ 

ff ^ = f dp' Ip' + |><P' - |l 


CO 


-^qP' 


(A.3) 

(A.4) 


where the first is in the position and the latter in the momentum basis. They have the properties 


A.+ A. Aj. A Iaa 

rj = r_^, Rl = R., (-Rx)' = l. (A.5) 

Their action on the position and momentum basis justifies their names 

Rxko) = 2 |2^ - ^o) RM = 2\2p- pf) (A.6) 


and 

f^\po) = \Po + fp) f^\qo) = \qo + (A.7) 

Moreover they form a group which is the representation of the affine group of reflections and 
translations, with the following composition laws 

44 = 4,4 = (A.8) 

tfR, = (A.9) 

They conform a pair of complementary orthonormal bases with the properties 

trfj7> = 2;rM(^-)8), trR^Ry = 2nhS{x-y), \xRfP^ . (A.IO) 

Switching to the double Dirac notation we rewrite the above as 

if^Wfp)) = 2nh6{^-fi), {{RMy)) = 2nhS(x-y), {(MU) = m-<’‘’f>.(A.ll) 

The labels x and f are related to the conjugate variables Q,P of ( [T7| ) as x = 2 and 
^ = JP. Thus we can think of reflection and translation operators as alternative position and 
momentum bases in double phases space. Using these properties we compute the quadruple 
traces needed in the main text: 

tr(ff,444) = 2nh6{,^, +^ 2 +^ 3 + (A.12) 

tr(Rxi^x,43^x4) = 2nndf^^^ - ^I^)^2<x,.x,>+2<x3,x4>_ (a.13) 
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We should remark at this point that traees of unitary operators are related semielassieally 
to elassieal periodie orbits and their aetions (21 0 , even when the elassieal and quantum 
evolutions are broken up into several steps [|47l . In the ease of four operators t^, sueh a 
trajectory is composed of four segments, each giving a phase space translation and forming 
a closed quadrilateral. The action of this trajectory is the symplectic area of the quadrilateral 
A 4 (^i,^ 2 ’^ 3 ’^ 4 ) = >) ^uits of h). In the case of reflections the 

trajectory connects the centers of the segments of this quadrilateral, which is a parallelogram 
Xi - X 2 + X 3 - X 4 = 0. Both cases are illustrated in Figure[^ 
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